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ABSTRACT 

Using the OSIRIS tunable narrow band imager on the 10.4m GTC (La Palma) we have 
mapped the SAB(rs)cd galaxy NGC 6946 over a ~7.3x7.5 arcminute 2 field in the emission lines 
of the [SII]AA 6717, 6731 doublet, and in Ha. From these maps we have produced catalogs of 
the Ha luminosities and effective radii of 557 HII regions across the disk, and derived the [SII] 
emission line ratios of 370 of these. The Ha observations were used to derive the mean luminosity- 
weighted electron densities for the regions of the sample, while the [SII] line ratios allowed us to 
derive values of the in situ electron densities in the denser zones from which the major fraction of 
the radiation in these lines is emitted, for 58 of the regions. This is by far the largest data set of 
its kind for a single galaxy. A classical two phase model is used to derive the filling factors of the 
regions. We find that although the mean electron density decreases with the square root of the 
radius of the regions, the in situ density is essentially independent of this radius. Thus the filling 
factor falls systematically, as the radius and the luminosity of the regions increases, with a power 
law of exponent -2.23 between filling factor and radius. These measurements should enhance the 
perspectives for more refined physical models of HII regions. 

Subject headings: galaxies: star formation — galaxies: individual (NGC 6946) — HII regions — galaxies: 
spiral 
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1. Introduction 

The initial theory of the ionization equilibrium 
in H II regions, proposed by Stromgrcn (1939), 
gave broad guidelines showing how their sizes and 
their emission line luminosities depend on the ion- 
izing luminosities of their central stars and the 
density of the surrounding ISM. It was not un- 
til twenty years later that Osterbrock and Flather 
(1959), measured the electron density in the Orion 
Nebula using the ratio of the [Oil] doublet lines 
AA3726,3729A and compared the result with a 
value of the rms electron density derived from the 
surface brightness in H/3, finding a difference of an 
order of magnitude. They inferred that this H II 
region must have a strongly inhomogeneous struc- 
ture. Assuming that the [Oil] line ratio is pro- 
duced in a distributed fraction of the gas having 
higher electron density (what we call the in situ 
electron density), they suggested that the frac- 
tional volume occupied of these denser zones, the 
" filling factor" could be found by assuming a sim- 
ple two phase model and comparing the two mea- 
sured densities. 

This basic scenario has been used in many sub- 
sequent studies of H II regions, both local and 
extragalactic. Simpson (1973) confirmed the de- 
ductions of Osterbrock and Flather (1959) for the 
Orion Nebula, and there have been many measure- 
ments of the in situ electron density in this, the 
nearest star forming zone, between then and the 
present day. We will cite here O'Dell & Goss who 
combined 327.5 MHz continuum observations with 
those of the optical [SII] doublet to explore spa- 
tial variations of the electron density and Rubin 
ct al. (2011) who used the [SII] doublet and the 
[SIII] infrared doublet from Spitzer spectroscopy, 
to derive values of the local electron density in this 
region of a few hundred per cm~ 3 . It is worth men- 
tioning a very different technique which yields the 
same qualitative result for HII regions. As an ex- 
ample Kassim et al. (1989), using only radio con- 
tinuum measurements at m and cm wavelengths 
also needed a highly inhomogeneous ISM to ex- 
plain their measured intensity ratios. More gener- 
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ally the rms electron density can be derived from 
the radio continuum, and a value of the in situ 
density from recombination line ratios in the opti- 
cal or near IR, with results for the latter which are 
canonically at least an order of magnitude greater 
than the former, pointing up the presence of the 
strong inhomogeneities in HII regions in general. 

The filling factor of the denser zones of ion- 
ized interstellar regions is of interest if we want 
to quantify the fraction of ionizing photons which 
escape from the environments of OB star clusters, 
either into the diffuse gas of the galaxy, or even 
into the intcrgalactic medium. It is well under- 
stood that an inhomogeneous hydrogen cloud is 
more permeable to Lyc photons. Following the in- 
ference by Ferguson et al. (1996) that escaping 
Lyc from H II regions in galaxies may be respon- 
sible for the ionization of the diffuse hydrogen well 
away from the regions, Zurita et al. (2002) pro- 
duced a model in which the surface brightness dis- 
tribution of the diffuse Ha in NGC 157 could be 
reproduced in detail in a model with specified Lyc 
escape factors from the regions. Giammanco et al. 
(2004) showed how inhomogeneities, if optically 
thick, can alter the derivations of all parameters: 
electron temperature and density and their vari- 
ances, ionization parameter, and abundance de- 
terminations. These considerations, not yet well 
enough studied, show the importance of quantify- 
ing density inhomogeneities. In consequence one 
would have expected to find in the literature many 
determinations of the key parameter: the filling 
factor. This is not, however, the case, basically be- 
cause the derivation of the in situ densities using 
suitable emission line ratios is not trivially easy. 
Detailed measurements with the the [SII] dou- 
blet AA6717,6731 emission line ratio using long slit 
spectroscopy were published for 9 luminous H II 
regions in 8 different galaxies by Castaheda et al 
(1992), but the most extensive data set is by Zarit- 
sky, Kennicutt & Huchra (1994) who measured 42 
individual H II regions in 39 galaxies. The range of 
those values, of order a few x 10 2 cm -3 , has been 
used repeatedly since then to make order of mag- 
nitude estimates of filling factors, by combining 
them with rms values of the electron density mea- 
sured using the radius and Ha surface brightness 
of individual HII regions, (as outlined below in sec- 
tion 3.2). Hunt & Hirashita (2009) presented [SII] 
derived in situ electron densities, based on HST 
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STIS data, for 24 HII regions, compact star form- 
ing zones, in a set of 24 galaxies, with correspond- 
ing rms densities obtained from H/3 luminosities 
and radii, and derived the filling factors for all of 
these. However the observations of Castaheda et 
al. (1992), of Zaritsky et al. (1994) and of Hunt & 
Hirashita (2009) do not permit a parametric study 
of behaviour within a single galaxy, to explore how 
the filling factor may vary systematically with e.g. 
the position of an HII region in the galaxy, or the 
overall physical parameters of the region itself: its 
Ha luminosity and its radius. 

In the present article we have used the capa- 
bilities of the OSIRIS tuneable narrow band filter 
on the 10.4m GTC (La Palma) to make narrow 
band images in the two emission lines of the [SII] 
doublet at 6716/6731 A across the full face of the 
SAB(rs)cd type galaxy NGC 6946, from which we 
can derive quite a wide range of in situ electron 
densities for individual HII regions with a range 
of luminosities and radii. At the same time we ob- 
tained flux calibrated images in Ha, from which 
we could derive values of the rms electron density. 
For those regions with reliably determined values 
of both of these parameters we could find their 
filling factors, and examine how those parameters 
vary according to the radius and galactocentric ra- 
dius of the regions. 

2. Observations and data processing 

The data employed in this study were obtained 
in an observing run during 2011 August 21 and 23, 
at the 10.4m Gran Telescopio Canarias (GTC), 
located at Roque de los Muchachos Observatory. 
The instrument used was OSIRIS, in tunable filter 
imaging mode, with a FOV of 7.3' x 7.5'. The mean 
seeing during the observations was about 0.7". We 
took a scan composed of 15 slices for Ha and 
the [Nil] AA6548, 6584 lines, tuned to wavelengths 
from 6549 A to 6633 A in the optical centre, with 
FWHMs from 13 A to 14 A. Each slice consisted 
of 6 exposures of 35 seconds, three on target and 
three off-target in order to sample the sky. For the 
[SII]AA6717,6731 lines, we took a scan composed 
of 19 slices, tuned to wavelengths from 6704 A to 
6812 A in the optical center, with FWHM~13A. 
Also, for these observations, each wavelength had 
6 exposures of 90 seconds, three on target and 
three off-target. In a forthcoming paper, we will 



discuss in detail the observation strategy and the 
reduction and calibration processes. As a short 
summary, the data were reduced and calibrated 
using the IRAF packageQ. Each image at each 
wavelength was bias-subtracted using bias frames 
taken during the night. The images for each wave- 
length were then flatfield-corrected using master 
flats obtained by the median combination of dome 
flatfields (five for each wavelength) and then cor- 
rected from ilumination by multiplying them by 
a normalized surface obtained from a twilight sky 
r' flatfield. Sky images for each wavelength were 
created by combining (with the median) the 3 off- 
target images after masking the stars. The sky 
image corresponding to a given wavelength was 
subtracted from the on-target images. After that, 
all the on-target images were aligned and an astro- 
metrical solution for each one was calculated. The 
position of the optical axis for each on-target im- 
age at each wavelength was determined using the 
positions of several image-ghost pairs. Employ- 
ing the program SExtractor (Bertin & Arnouts, 
1996) and the sum of all Ha wavelengths as the 
detection image, a catalogue of H II regions for 
each wavelength was created. A region was con- 
sidered detected when it presented a flux per pixel 
3 sigma over the background, and when the total 
area (defined as the isophotal area) was larger than 
the area subtended by a circle with radius equal 
to the seeing of the image. The equivalent radius, 
r, for each H II region is defined as the square root 
of the isophotal area of the region divided by 7r. 

For each H II region detected in each tuned 
slice, the corresponding real wavelength A was cal- 
culated, following Gonzalez et al. (2013, in prepa- 
ration^]), as 

A = A - 5.047-2 + a 3 (A)ri (1) 

Where Ao is the central wavelength of a given slice, 
r a is the distance of the H II region to the optical 
centre in arcmin, and 03(A) is defined as 

a 3 = 6.17808 - 1.6024 x 10~ 3 A + 1.0215 x 10~ 7 A 2 

(2) 



X IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Uni- 
versities for Research in Astronomy (AURA) under coop- 
erative agreement with the National Science Foundation. 

2 For more information, see the webpage 

http://gtc-osiris.blogspot.com.es/2012/03/new-rtf-calibration.html 
and follow the links there. 
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For each region detected, the pseudospectra of 
the Ha zone and the [SII] zone were generated. 
The data were then flux calibrated following the 
prescription by Jones et al. (2002) for an equiva- 
lent system to that used here. 

In order to obtain the amplitudes of the ex- 
pected components, we performed a least squares 
inversion of the spectral profile. The method uses 
a combination of Gaussians whose number and 
widths are pre-assigned according to the expected 
emission lines, leaving their amplitude and cen- 
tral wavelengths to vary freely inside a range. The 
Gaussians are then convolved with the spectral re- 
sponse function of the instrument (an Airy func- 
tion, see Jones et al. 2002); solving the least- 
squares problem yields the amplitudes and central 
wavelengths of all the components. 

The errors for each line were calculated as the 
addition of the uncertainty in the measurement 
of the total integrated flux of the region and the 
uncertainty in the calibration. The error in the fit- 
ting process is small enough, when compared with 
the previous ones, to be considered negligible for 
the calculation. The mean errors in the line inten- 
sities are 5% for Ha and 7% for the [SII] doublet. 

In Fig. [T]we show the resultant image of adding 
all the Ha wavelengths (upper panel). A typ- 
ical H II region with enough signal to noise in 
Ha and [SII] (RA: 308.6998566 degrees, DEC: 
+60.1772541 degrees) was selected as an example 
and its position is indicated by a circle in the im- 
age. In the lower panels is represented the spectral 
energy distribution of the selected region in the 
Ha+[NII] zone (left panel) and in the [SII] zone 
(right panel). The emission in each wavelength is 
indicated by the open circles, the dashed lines are 
the best fit to the data, and the continuous lines 
are the deconvolved components for each emission 
line. 

The final catalogue has 557 H II regions with 
Ha emission detected and measured, with 370 of 
these H II regions having measured fluxes in both 
of the [SII]AA6717,6731 lines. In Fig. His repre- 
sented the numbers of detected regions in bins of 
0.1 R/R25 (with i?25=5.7', from de Vaucouleurs et 
al. 1991) for Ha (empty bars) and for [SII] (filled 
bars). 



3. In situ and mean electron densities 
3.1. In situ electron density 

The lines of R [S/7] = [SII] A6717/[SII]A6731 ratio 
are emitted by different levels with similar exci- 
tation energy, so this ratio is determined by the 
collisional strengths, and it is, therefore, sensitive 
to the electron density (Osterbrock, 1989). To ob- 
tain the electron density from R^//], we used the 
analytical solution given by McCall (1984) for a 
three-level atom: 



R 



[SII] 



1.49 



1 + 3. 77 a; 



1 + 12.8 a; 



(3) 



where x = 10~ 4 n e i~ 1/2 , t = T/10000, n e is the 
electron density in cm -3 , and T is the electron 
temperature in K. If we assume that t — 1 for all 
the regions, this gives us: 



1.49 R 



[SII] 



12. 8R 



[SII] 



5.6713 



x 10 4 



(4) 



According to Castaheda et al. (1992) or Zarit- 
sky et al. (1994), the lower limit for reliably mea- 
sured densities is 10 cm' 3 . However, due to the 
uncertainties in the fluxes for the [SII] lines, and 
the functional dependence of the density on the 
line ratio, we have selected a more conservative 
value of 50 cm~ 3 as a lower limit for the reliably 
measured density. With this limit, we were able 
to obtaint the electron density for 58 regions, well 
distributed within the disk of the galaxy. 

In Fig. [3]we have represented the electron den- 
sity as a function of the equivalent radius for the 
H II regions. The values of the density range 
from ~50cm -3 (our imposed lower limit) up to 
~10000cm -3 , with a median value of 279cm -3 . 
There is no clear correlation betweeen these two 
parameters. 

3.2. Mean electron density 

The mean luminosity- weighted electron density 
< n e >, is determined from the Ha luminosity 
and the equivalent radius r of each H II region. 
Following Gutierrez & Beckman (2010) (there- 
after GB10), this mean electron density can be 
expressed as: 



< n e >= 1.5 x 10 



L(Ha) 



(5) 
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Where L(Ha) is the luminosity in the Ha line and 
r is the equivalent radius of the H II region. If 
L(Ha) is expressed in ergs -1 and r in pc, eq. [5] 
gives < n e > in cm -3 . The values found for < 
n e > lie in the range a few cm -3 to a few tens 
cm -3 (see also Rozas et al., 1996, GB10). 

In Fig. 0] we represent the logarithm of the 
mean electron density as a function of the log- 
arithm of the equivalent radius for all the H II 
regions with Ha data (left panel) and for the re- 
gions where we were able to determine the in situ 
electron density. We have corrected the mean elec- 
tron density for the effect of the dependence on the 
galactocentric radius using the method proposed 
in GB10. A least squares fit to the data, gives the 
result 



< n P >= 



51 



,0.55±0.05 ' 



(6) 



which implies that the larger regions have lower 
mean electron densities. This result agrees with 
that presented in GB10, where they find that the 
relationship between the equivalent radius of the 
H II regions and their mean electron density in 
two quite different galaxies: M51 and NGC 4449, 
follows the expression 



1 

r 0.5 ' 



(7) 



4. Filling factor 



We assume for an H II region that the emis- 
sion line ratios detected are weighted by density, 
as they are formed principally within the denser 
zones which occupy just a fraction of the total vol- 
ume of the region, and that the rest of the H II re- 
gion is formed by a low density component. This 
component is sampled, at the same time as the 
high density clumps, when we derive mean elec- 
tron density (Rozas et al. 1996, GB10). Following 
Osterbrock & Flather (1959), we define the filling- 
factor (//) of an H II region as 



< n e . > 



(8) 



In Fig. [5] we present the logarithm of // as a 
function of the logarithm of the equivalent radius 
of the H II regions of NGC 6946. A clear anticorre- 
lation exists between these parameters: larger H II 



regions have lower values of the //. The correla- 
tion coefficent for this relationship is -0.3822, with 
a value of p=0.0004, where p is the probability of 
getting a correlation as large as the observed value 
by random chance. In this case, p is small enough 
(about 0.04% ) to conclude that the two quantities 
are related. A first degree polynomial (shown in 
the figure as a continuous line) was fitted to the 
data using the least squares method, and can be 
expressed as 

This relationship indicates that, while the dense 
clumps within an H II region maintain a density 
independent of the luminosity of the region (see 
Fig. [3]), for larger H II regions, there is a larger 
volume of low density gas filling the space between 
the denser zones. 

5. Conclusions 

We have obtained reliable luminosities and radii 
for 557 H II regions of NGC 6946 in Ha and 370 of 
these H II regions in both of the [SII] AA6717,6731 
doublet lines. 

From the doublet line ratio we could measure 
reliably in situ electron densities, n e , for 59 H II 
regions. We found no correlation between the in 
situ electron density and the equivalent radius of 
the H II regions. This seems to indicate that the 
denser zones where the bulk of the [SII] line emis- 
sion is produced, show little variation in density 
over the full range of H II region sizes and lumi- 
nosities. 

We also obtained the mean electron density < 
n e > for all the regions, using Ha luminosities and 
the effective radii. We found, as also found by 
GB10, a power law relationship with an exponent 
of -0.55, between < n e > and the equivalent radius 
of the H II regions. 

We then calculated the filling factor for the H II 
regions of the sample. We found a power law rela- 
tionship (with an exponent of -1.78) between the 
filling factor and the equivalent radius of the H II 
regions, which indicates that larger H II region 
are more porous. We can explain our results with 
the scenario in which the denser zones, with elec- 
tron densities in the range 100 - 1000cm -3 , do 
not change their internal properties as we go from 
smaller, less luminous to large, more luminous 
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regions, but they occupy a progressively smaller 
fraction of the volume. The tenous interclump gas, 
with density in the range 1-10 cm~ 3 , occupies an 
increasing fraction of the volume the larger and 
more luminous the H II regions are. We have pro- 
duced the first quantitative estimate of this be- 
haviour, which should open the way to more re- 
fined physical models of H II region structure. 
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Fig. 1. — Upper panel: composite image of NGC 6946 obtained by adding all the wavelengths for the Ha 
scan. The circle indicates the position of region at RA: 308.6998566 degrees, DEC: +60.1772541 degrees. 
Lower panel: spectral energy distribution for Ha + [Nil] (left) and [SII] (right) zones. Open circles are the 
fluxes measured in each wavelength, the dashed line is the best fit for the data points, and the continuous 
lines are the deconvolved, deblended, and continuum subtracted gaussians for each emission line. 
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Fig. 2. — Number of detected regions as a func- 
tion of the galactocentric radius, in bins of 0.1 in 
i?/i?25- The numbers of regions detected in Ha 
are the empty bars, the filled bars are the regions 
detected in [SII] . 
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Fig. 3. — Logarithm of the in situ electron den- 
sity as a function of the equivalent radius for the 
measured H II regions. 
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Fig. 4. — Logarithm of the mean electron density, corrected by a factor which takes into account the 
galactocentric distance, as a function of the logarithm of equivalent radius. The continuous line is the best 
fit for the data points. The dashed lines indicate the 2a confidence for the fit. In the left panel all the H II 
regions with Ha data are represented. In the right panel, only the regions with a determination of the in 
situ electron density are represented. 
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Fig. 5. — Logarithm of the filling factor as a func- 
tion of the logarithm of the equivalent radius of 
the regions. The continuous line is the best fit for 
the data points. The dashed lines indicate 66% 
confidence level for the fit. 



